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1 INTRODUCTION 



A detailed study of the ionization history of the Universe 
is fundamentally important for our understanding of the 
properties of the structure and evolution of the Universe, 
particularly the large-scale structure and galaxy formation. 
Although the epoch of galaxy formation is often referred as 
the dark age due to the difficulty in direct observations, it 
is, nevertheless, feasible to investigate in details the ioniza- 
tion history of the Universe through the cosmic microwave 
background (CMB) anisotropics and polarization. The re- 
cent CMB experiments, such as the BOOMERANG (de 
Bernardis et al. 2000), MAXIM A-1 (Hanany et al. 2000), 
CBI (Mason et al. 2002), VSA (Watson et al. 2002), DASI 
(Halverson et al. 2002), and its polarization data (Kovac 
et al. 2002; Leitch et al. 2002) have shed light on probing 
the dark age of the Universe. In particular, the newly re- 
leased WMAP temperature-polarization data (Hinshaw et 
al. 2003; Kogut et al. 2003) provide a new way for under- 
standing of the very early stages of galaxy and star for- 
mation. We expect to have the future Planck polarization 
data with unprecedented accuracy. The polarization power 
spectrum from these two missions will therefore provide us 
the information about the kinetics of hydrogen recombina- 
tion and allow us to determine the parameters of the last 



ABSTRACT 

We investigate two sets of two-epoched reionization models and their manifestation 
in the CMB anisotropy and polarization data of the recent WMAP project and make 
some predictions for the future Planck missions. In the first set of models, the uni- 
verse was reionized twice, first at z ~ 15 by population III stars and then at z ~ 6 by 
stars in large galaxies. In the second set of models, the extra peak-like reionization at 
high redshifts z > 100 is induced by the decay of unstable particles, followed by the 
standard picture of reionization at z ~ 6. We examine the general properties of these 
two-epoched reionization models and their implication in the WMAP CMB tempera- 
ture anisotropy-polarization cross-correlation. We have shown that these models have 
comparable likelihood values for the WMAP data and distinct characters which can 
be tested with the Planck high sensitivities. 
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scattering surface and the ionization history of the cosmic 
plasma at very high redshifts z ~ 10^. 

In the framework of the modern theory of the primary 
CMB anisotropy and polarization formation, the theory of 
hydrogen recombination are assumed to be a 'standard' 
one. The classical theory of hydrogen recombination for the 
pure baryonic cosmological model was developed by Pee- 
bles (1968), Zel'dovich, Kurt and Sunyaev (1968), and was 
generalized for non-baryonic dark matter by Zabotin and 
Naselsky (1985), Jones and Wyse (1985), Seager, Sasselov 
and Scott (2000), Peebles, Seager and Hu (2000). This stan- 
dard model of recombination has been modified in various 
ways. 

First of all, there are some variants from the standard 
hydrogen recombination model, namely, the delay and ac- 



celeration of recombination at the redshift Zr, 



lO"* due to 



energy injection from unstable massive particles (Doroshke- 
vich and Naselsky 2002) or due to the lumpy structure of 
the baryonic fraction of the matter at small scales (Naselsky 
and Novikov 2002) , in which the typical mass of the clouds 
is of the order lO'^ - 10*^ Mq (see Doroshkevich et al.2003 
and the references therein). Secondly, the most crucial part 
of the ionization history of the Universe is related to the 
large-scale structure and galaxy formation and is called late 
reionization. The model of the late reionization is not yet 
well-established and needs further investigations. 
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The conventional view of the ionization history* is that 
cosmological hydrogen became neutral after recombination 
at Zrec — 10^ and was reionized at some redshift Zreion, 
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where Tr is the Thomson optical depth, Qh is the present 
baryonic density scaled to the critical density, 0dm is the 
dark matter density, h = J7o/100kms~'^Mpc is the Hubble 
constant, (Yp) is the helium mass fraction of matter. Re- 
cently Cen (2002) has proposed the model of the late reion- 
ization with two epochs. Firstly, hydrogen was reionized at 
redshift zUl^^ ~ 15 by Population III stars and secondly at 

(2) 

•^reion — 6 by stars in large galaxies. On the other hand, we 



also discuss another distinct feature of reionization model, 
which is called the peak-like reionization. This shoot-up in 
the ionization fraction at z > 100 can be induced by energy 
injection into the cosmic plasma. 

These two-epochcd reionization models, which can be 
tested by the WMAP and future Planck data (Cen 2002), 
would be significant for the interpretation of the polarization 
measurements. The polarization of the CMB from the late 
reionization epoch (or epochs) is sensitive to the width of 
the period A^reion, when the ionization fraction Xe increases 
from the residual ionization {xe ~ 10~^) up to ~ 0.1 — 1 
(Seljak and Zaldarriaga 1996). They can provide unique in- 
formation about the physical processes induced by compli- 
cated ionization regimes. The aim of the paper is to discuss 
the distinct characters of these models in the light of recent 
WMAP data and to predict the peculiarities in the polar- 
ization power spectrum induced from both the Cen model 
(2002) of the late reionization and the extra peak-like reion- 
ization, taking into account the properties and the sensitiv- 
ities of the upcoming polarization measurements. 



2 PHENOMENOLOGY OF THE 

TWO-EPOCHED LATE REIONIZATION 

The model of the reionization process proposed by Cen 
(2002) can be described phenomenologically in terms of the 
injection of additional Ly-c photons via the approach by 
Peebles, Scager and Hu (2000), Doroshkevich and Naselsky 
(2002), Doroshkevich ct al.(2003). For the epochs of reion- 
ization the rate of ionized photon production m is defined 
as 



drii 
IF 



= Si{z)nb{z)H{z), 



(2) 



where H{z) and nh{z) are the Hubble parameter and the 
mean baryonic density at z, respectively, ei(z) is the effec- 
tiveness of the Ly-c photon production. As one can see from 
Eq.(??) the dependence of ei{z) parameter upon redshift z 
allows us to model any kind of ionization regimes, including 



In this paper we refer to "reionization epoch(s)" as epoch(s) 



with Xe above the residual fraction (^ 
the recombination epoch. 



10 ) of ionization from 



heavy particle decays. This parameter also includes uncer- 
tainties from the fraction of baryons that collapse and form 
stars, and the escape fraction for ionizing photons. For late 
reionization, the ionization fraction of matter Xe = Ue/n 
can be obtained from the balance between the recombina- 
tion and the ionization process 



(3) 



s cm ' is the 



^ = -areciT)nhx'i + ei(2)(l - Xe)H{z), 

where Or^dT) ~ 4 x lO"^-"* (r/10*A') 
recombination coefficient and T is the temperature of the 
plasma and rib is the mean value of the baryonic number 
density of matter. In an equilibrium between the recombi- 
nation and the ionization process the ionization fraction of 
the matter follows the well-known regime 



xljz) 
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where H{z) = i?oi/Oni(l + z)^ + l — and Ub — 2 x 
10"''(nb/i^/0.02)(l -I- z)^. We would like to point out that 
Eq.(4) can be used for any models of the late reionization, 
including the Cen model (2002) by choosing the correspond- 
ing dependence of the ei{z) parameter on redshift. This 
point is vital in our modification of the RECFAST and the 
CMBFAST packages, from which we can use the standard re- 
lation for matter temperature T{z) ~ 270 (1 + z/lOO)^ K 
and all the temperature peculiarities of the reionization 
and clumping would be related with the ei(z) parameter 
through the mimic of ionization history. For example, in the 
Cen model (2002) the function T{t) has a point of maxima 



J- max 

a,t z < z. 
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1.5) X 10'* at 2 ~ z^lion ^^'^ decreases slowly 
down to T{t) ~ 10* ~ const at the redshift 
range 6 < « < 12. Let us introduce some model of the ei{z) 
parameter dependence over z as 



ei{z) = So exp 



+eiil+z)-"'Q{z 



(1) 



,(5) 



where eo ,£i and m are the free parameters, Azi ^ ^reion 
is the width of the first epoch of reionization and Q{x) is 
the step function. The first term of Eq.(5) corresponds to 
peak-like reionization at 2 = Zreion, which decreases signifi- 
cantly at z > ^reion- The second term is related to modelling 
the second epoch of the reionization model discussed by Cen 
(2002), which results in a monotonic increasing in ei{z) func- 
tion as a function of time. From Eq.(5) at a ~ ^^eion ^® 
obtain 



Qlrec ( ^reion ) '^b ( -^reion ) 



(1) 
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where 



eo > \zllljarcc{zilljnb{zlllj 
u2 \ -1/2 
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One can find from Eq.(4) and (6) that for the second epoch 



of reionization {xe ~ 1 at 2; ~ 



(2) 



6 in the Cen model 



(2002)) the amplitude of the ei parameter and the index m 

should satisfy the following 



(8) 
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Figure 1. The ionization fraction for the two-epoched reioniza- 
tion models. Solid line corresponds to the model 1, the dash line 
the model 2 and the dash-dot line the model 3. 

where 7 is the parameter between 0.1 to 1 in order to model 
the properties of the Ccii model and its variants. We demon- 
strate the effectiveness of our phenomenological approach in 
Fig. 1: the ionization fraction Xe against redshift for the 
three models listed below: 

• model 1: eo = 1.3 x 10^ ei = /3 x 10^ /3 = 1, m = 7; 

• model 2: eo = 1.3 x 10^ £1 = ^ x lO'', /? = 0.1, m = 7; 

• model 3: eo = 1.3 x 10^el =px 10^, /3 = 100, m = 8. 

The curves are produced from the modification of the rec- 
FAST code (Seager, Sasselov and Scott 2000). Obviously the 
model 2 (dash line) is not physically viable, we neverthe- 
less can use it to test the sensitivity of CMB polarization to 
non-monotonic shape of ionization fraction. For all models 
we use the following values of the cosmological parameters: 
fibft^ = 0.022, n^h"^ = 0.125, nx = 0.7, h = 0.7, n^ + D.x = 
1. We show the 3 models in Figure 1. Model 1 and 2 mimic 
the properties of the Cen model (2002) where there is a dip 
in the reionization fraction, whereas model 3 corresponds to 
roughly the standard reionization model with no significant 
presence of peaks in ionization fraction. 



3 THE CMB POLARIZATION FOR THE 
TWO-EPOCHED LATE REIONIZATION 

MODELS 

In order to find out how sensitive the polarization power 
spectrum is to the two-epoched reionization models, we con- 
sider phenomenologically the different variants of hydrogen 
reionization models by modifying the RECFAST (Scagcr, Sas- 
selov and Scott 2000) and CMBFAST code (Seljak and Zaldar- 
riaga 1996). In Fig. 3 we plot the polarization power spec- 
trum Cp{(.) for the model 1, 2, the standard single reion- 
ization model at Zrdon — 6 and 13.6. The difference be- 
tween model 1 and 2 mainly lies in between the multipoles 
2<£<30. 

3.1 The anisotropy and polarization in 
comparison with the WMAP data 

To characterize the differences between the reionization 
models we have used the WMAP programme (Verde et 



model variants 12 3 



Likelihood (T) -494.357 


-525.928 


-495.159 


Likelihood (TE) -231.601 


-223.749 


-230.766 



Table 1. The likelihood parameters of the variants of the Cen 
model (2003). 
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Figure 2. The anisotropy power spectrum for models 1 — 3 in 
comparison with WMAP and CBI data. All the models predict 
approximately the same shape which cannot be seen in the fig- 
ure due to degeneracy. The corresponding likelihood parameters 
are shown in Table 1. The WMAP data are taken from the offi- 
cial website at the range £ < 500. The CBI data are taken from 
CBIMl and M2 data sets. 



al. 2003) for the calculations of the likelihood for the 
anisotropy and TE correlation power spectra against those 
from WMAP results, shown in Table 1. As one can see, for 
all the models of late reionization we get fairly good con- 
sistency for anisotropy and the TE correlations. The accu- 
racy of the WMAP data (Hinshaw et al. 2003; Kogut et 
al. 2003) however is not enough for discrimination between 
the models. Therefore, the future Planck data is required 
for more accurate investigations of the history of hydrogen 
reionization at relatively low redshifts z < 30. For compari- 
son, the WMAP best-fitting cosmological model has the like- 
lihood parameter —486.245 (with 900 data points) for the 
anisotropy and —228.695 (with 499 data points) for the TE 
cross-correlation. As one can see from Table 1, all the mod- 
els have excellent agreement with the WMAP data, but we 
are not able to distinguish with the accuracy of the WMAP 
data between the two-epoched reionization models in the 
Con model (2002). In order to check the high multipole 
range of the power spectrum of temperature anisotropics 
in the models 1 — 3, in Fig. 2 we plot the C{1), which are al- 
most indistinguishable for all the models and are consistent 
with the WMAP and the CB/data. The result shows excel- 
lent agreements between the theoretical curves and the data 
points from the experiments. The most intriguing question 
is could the future Planck data sets allow us to distinguish 
any peculiarities of the late reionization epochs? 
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Figure 3. The polarization power spectrum for different models 
of the reionized universe. The solid line corresponds to the model 
1, the dotted line the model 2, the dash and the dash-dot line are 
the model with single reionization (a^eCzreion) = 1) at ^reion — 6, 
and at Zreion — 13.6, respectively. 



Figure 5. The plot of -D^j, the deviation in polarization for 
different reionization models. The solid line D^,! corresponds to 
the models 1 and 2, the dash-dot line D^ i is from the model 1 in 
comparison with the model of single reionization at jSreion — 

13.6 

and the dash line Ds,i is from the model 1 in comparison with 
the model of single reionization at Zroion — 6. The thick solid 
(binning) and dotted curves (/sky = 0.65, no binning) mark the 
estimated Planck errors. 
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Figure 4. The TE cross-correlation power spectrum for different 
models of the reionized universe. The solid line corresponds to 
the model 1, the dotted line the model 2, the dash and the dash- 
dot line arc tlie model with single reionization (xe(2^reion 

) = 1) at 

^reion — 6, and at Zreion — 13.6, respectively. 



3.2 The anisotropy and polarization in 

compcirison with the future Planck sensitivity 

The differences between the late reionization models in com- 
parison with the expected sensitivity of the Planck mission 
can be expressed in terms of the power spectrum Cp{tj (for 
the anisotropy, E and TE component of polarization) 

where the indices i and j denote the difFcrcnt models. 

In order to clarify the manifestations of the complex 
ionization regimes in the models 1 and 2 we need to com- 
pare the peak to peak amplitudes of the Dij{i) function 
with the expected error of the anisotropy power spectrum 
for the Planck experiment. We assume that the systematics 
and foreground effects are successfully removed. The corre- 
sponding error bar should be 

TOT! 



where w = (crp6lFWHM) ^, Wi ~ exp [-£{£ + l)/2£i] , /sky ~ 
0.65 is the sky coverage during the first year of observations, 
(Tp is the sensitivity per resolution element Sfwhm x ^fwhm 

and is = \/81n2ep.^jjM- 

For all Planck frequency channels, for example, the 
FWHM are less than 30 arcmin, so for the estimation of 
the errors at £ < 40 range we can omit the second term in 
Eq. (10). In Fig.3 we show the polarization power spectrum 
and TE cross-correlation for models . While polarization 
power spectrum is not observed by WMAP, the upcoming 
Planck mission will be able to provide us the observation to 
differentiate different models. It is clear from Fig. 4 that the 
TE cross-correlation spectra for different regimes of reion- 
ization have different shape at £ < 10 range but do not vary 
significantly for £ > 10. All the deviations lie inside the cos- 
mic variance and are practically not observable. 

As one can see from Fig. 5 for 73i,2(^) the corresponding 
peak to peak amplitudes are on the order of magnitude 20% 
at ^ ~ 10-40, while the errors AC{£)/C(£) are in about the 
same one. Such small deviations in the polarization power 
spectrum caused by the complicated ionization regimes can 
not be tested directly for each multipole of the C(Z) power 
spectrum by the Planck mission, even the systematic effects 
would be removed down to the cosmic variance level. 

As shown in Fig. 5, the deviation D2,i mostly lies inside 
the error region. This indicates that the upcoming Planck 
observational data would not be able to distinguish the two- 
epoched late reionization models from each other, where the 
only difference is in the amplitudes of the minima of ion- 
ization fraction. However, it is worth noting that both the 
models 1 and 2 have significant deviation from the standard 
single reionization model (the dash and the dash-dot line). 

The shape of the polarization power spectrum in the 
two-epoched reionization model differs from the shapes for 
single reionization models even for ^reion — 13.6. Such kind 
of dependence is related to the difference with Az of the 
epoch when the ionization fraction starts to grow from 

Xe,min ~ 10 Up tO Xg^max ^ 1. 

An unique possibility to detect more complicated struc- 
ture of late reionization would be from the binning of the ini- 
tial data, using, for example, the same range of bin ~ 15—30. 
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In such case, if the correlations between each multipole are 
small, the accuracy of the C{£) estimation would be approx- 
imately 4 — 5 higher then for unbinning power spectrum and 
non-monotonic structure of the D{£) function could be de- 
tectable for the anisotropy and for E and TE polarization 
as well. 



4 PEAK-LIKE REIONIZATION AT HIGH 
REDSHIFTS 

In this section we shall investigate another type of two- 
epoched reionization models. What is the implication of an- 
other reionization occurring at high redshifts, if, for exam- 
ple, one of the epochs of the pro-roionization took place at 
redshifts 30 < 2 < 1000? Note that for z > 30 the Comp- 
ton cooling of the plasma is extremely important and any 
energy injection to the cosmic plasma could produce rela- 
tively short epochs of reionization, when the ionization frac- 
tion became significantly higher (xe,max ~ 1), but for rel- 
atively short time interval. We call such distinct character 
of reionization the peak-like reionization. Such regimes can 
be induced by the decay or annihilation of some unknown 
particles ore decay of the primordial black holes (Naselsky 
1978; Ivauov, Naselsky and Novikov 1994; Kotok and Nasel- 
sky 1998) during the long period 3 x 10^ — 10** years. Be- 
cause of the Compton cooling of the plasma the injected 
energy density einj would be absorbed by the CMB photons 
leading to j/-distortion in the black-body spectra. Peebles, 
Seager and Hu (2000) have shown that the corresponding 
value of the y-parameter in this model is y ^ 0.25einj/ecMB, 
where ecMB is the energy density of the CMB at the red- 
shift of the injection. Taking into account the COBE upper 
limit for y-parameter (Fixen et al. 1996) Xjcobe < 2 x 10~^, 
one can estimate the upper limit of the energy injection 
£inj < ^Vcobe^CMS- On the other hand, for reionization of 
each hydrogen atom we need to have roughly one photon 
with energy E ~ I, where I ~ 13.6eV. Thus, einj ~ Xelub 
and we obtain the limit 



/ ecMB 



^43/c, 



10 (l + -inj)(Q^ 



£CMB 
2 



(11) 



where rup ~ 1 GeV is the proton mass and eb is the energy 
density of baryons at rcdshift Zinj- From Eq.(ll) it is clear 
that the peak-like reionization [xe ~ 1) is self-consistent 
with the COBE observational limit on j/-parameters. We can 
describe the peak-like reionization in terms of the injection 
of an additional Ly-c photons as in Section 1, but with the 
source term now the Dirac (5-function 



^ =r]nb{z)5D{t-tp), 



(12) 



where q is the effectiveness of the Ly-c photon production, 
and tp is the age of the Universe at the moment of peak-like 
reionization. Thus for the ionization fraction Xe = we get 



,{T)nbxt + vC^ - Xe)5{t - tp). 



(13) 



dXe _ 
It " ~ 

Quantitatively we can assume that for the first reioniza- 
tion epoch the maximum of the ionization fraction is small 
{xe 1) and we neglect Xe in the 1 — Xe term in Eq.(13). In 
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Figure 6. The ionization fraction for the peaJi-like reionization 
model (model 4). The dotted, the solid and the dash-dot lines 
correspond to model 4a, 4b and 4c. 



such an approximation the ionization balance of plasma is 
determined by the pure balance between the recombination 
term and the energy injection term in Eq.(12): 



Xe{t) ~ 7? 1 -h?? / a{T)nbdt 



\ 



(14) 



/ 



for the temperature of the plasma we can assume that 
T{tp) ~ (1 - 2) X lO"* K and at t > tp, 



dT _ TcMB — T 
dt Tc 



(15) 



where Tc is the Compton cooling time. Note that our as- 
sumption about maxima of T{tp) ~ 10~^K is closely related 
with the energy balance. For reionization of the hydrogen 
by electromagnetic cascades the typical energy release is in 
the order luhXe. Such a part of the energy of a cascade does 
not behave as a simple (5-function in the energy spectrum 
(see Doroshkevich & Nciselsky 2002) and is usually charac- 
terised by power law E ' . Because of Compton scattering 
of high energy 7-quanta off electrons, it is natural to esti- 
mate the kinetic energy of electrons as UhkT ~ In^Xe and 
T ~ Ixe/k ~ lO'^K whore k is the Boltzmann constant. At 
z > 4Q the ratio between characteristic time of recombina^ 
tion tree = [ot{T = 10''K)nb]~^ and Tc is more than one 
order of magnitude, while both of them are much smaller 
than the Hubble time. The relaxation of the matter tem- 
perature to the CMB temperature proceeds faster than the 
ionized hydrogen becoming neutral. Thus, while the temper- 
ature of matter is close to the CMB temperature Tcmb , the 
corresponding time of recombination is 



Mr. 



\dXe/dt\ 



; ?7 ^tr(T'cMB), 



(16) 



where tr = tree at T = Tcmb. Taking into account the 
above-mentioned properties of the temperature history of 
the plasma we can estimate the Thomson optical depth Ar^ 
caused by the peak-like reionization as follows 



In 



(17) 
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Figure 7. The polcirization power spectrum for the extra peaJc- 
like reionization models (the model 4). The dash, the solid and 

the dash-dot line correspond to the model 4a, 4b and 4c, respec- 
tively. The dotted line is the standard single reionization model 

at S^reion — ^• 



Figure 8. The polarization power spectrum for the model 4 as in 
Fig. 7 but in logarithmic scale, in which the differences between 
the power spectra can be seen more clearly at the multipole range 
t < 100. 



where Tr satisfies Eq.(l) at «reion = «reion(tp). For example, if 
the peak-hke reionization takes place at the redshift Xreion = 
200 and 77 ~ 0.1 - 0.2 then 



^ -0.1(1 + 2,, 



0.022 



0.125 



-1/2 



(18) 



and the corresponding values of the optical depth from 
Eq.(17) are in order of the magnitude Ar^ ~ 0.03 — 0.05. As 
one can see, if the CMB anisotropy data are consistent with 
the limit on Thomson optical depth of reionization r,- < 0.1 
(Doroshkovich, Naselsky, Naselsky and Novikov 2003), then 
roughly 30 — 50% of Tr can be induced by the peak- like reion- 
ization and 50 — 70% with the late reionization caused by the 
structure formation at low redshifts. To describe the peak- 
like reionization numerically we use a Gaussian approximar 
tion for the energy injection in Eq.(12) 



ri6{t — tp) — » £,H{z) exp 



(A2) 



(19) 



and describe the following two sets of examples for the peak- 
like reionization: 



model 4a: 2reion ■ 
model 4b: 2reion 

model 4c: Zreion '■ 



and 



model 5a: 2reion ■ 
model 5b: ^reion ■ 
model 5c: 



: 200 with : 
: 200 with Az : 
200 with A2 -- 



: 500 with Aa : 
: 500 with A2 : 
500 with A2 = 



5, ^ = 100; 
= 5, 5 = 10; 
5,^ = 1, 



: 12.5, C = 100; 
: 12.5, 5 = 10; 
12.5, 5 = 1, 



where Az/2reion = 0.025 for both sets of models. The late 
reionization a::e(z ~ 6) = 1 at 2 ~ 6 is included in both 
model 4 and model 5, as it is the standard part of the CMB- 
FAST package. In Fig. 6 wc plot the shape of the ionization 
fraction of the plasma for the model 4. As an analyti- 
cal description, one can see the peaks of Xe at 2reion = 200, 
which drops down at 2; ~ 100 — 150. 

By modification of the CMBFAST code for the primary 
polarization, we plot in Fig. 7-9 the corresponding power 
spectrum for the model 4 and 5. As one can see from Fig. 7 
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Figure 9. The polarization power spectrum for the extra peak- 
like reionization models (the model 5) in logarithmic scale. The 
dash, the solid and the dash-dot line correspond to the model 
5a, 5b and 5c, respectively. The dotted line is the standard single 



reionization model at Zr, 



■ 6. 



the common manifestation of the extra peak-like and the 
standard late reionization produce interesting features in the 
power spectrum. Namely, for high multipoles i the ampli- 
tude of the power spectrum decrease as exp(— Tr), while for 
the multipole range I < 100 the manifestation of the peak- 
like reionization is very clear, in particular at 10 < ^ < 200 
from Fig. 8. Figure 9 shows how important the peak-like 
reionization at Zrcion = 500 could be for the distortion at 
higher multipole range of the polarization power spectrum. 
In Fig. 10 we plot the D^^a function, the comparison between 
the model 4 and the standard single reionization model at 
2:rcion — 6. Fig. 1 1 is thc comparison between the model 5 
and thc standard single reionization model at Zrcion — 6. In 
Tabic 2, wc show the calculated likelihood for the anisotropy 
and TE correlation power spectra from the model 5 against 
those from WMAP results (Hinshaw et al. 2003; Kogut et 
al. 2003). They are close to the parameters from the model 
4. 

Once again we would like to point out that all the pe- 
culiarities induced by the extra peak-like reionization have 
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Filter peak-like model 2 



model variants 5a 5b 5c 



Likelihood (T) 


ruled out 


-770.693 


-509.583 


Likelihood (TE) 


-341.908 


-261.204 


-236.388 



Table 2. The likelihood of the Vciriants of the peak- like model 
(model 5). 



localized structure which appears at some fixed multipole 
range. These features can by tested by the Planck polariza- 
tion measurements. 



5 CONCLUSIONS 

We have investigated the two-epoched reionization models 
of the Universe. The two-epoched reionization can be in- 
duced by the structure formation as described in Cen model 
(2002), or caused by unknown sources of the energy injection 
(peak- like reionization) at relatively high redshifts z > 30. 
We have shown that for the Cen model (2002) the WMAP 
and the Planck mission would be able to detect the general 
shape of the ionization history for the two-epoched reionized 
plasma, which differs from the single reionization models at 
z ~ 13.6 or z ~ 6. However, any peculiarities of the ioniza- 
tion fraction of the matter inside the range 6 < z < 13.6, 
such as the decreasing of ionization, do not observed by the 
WMAP experiment due to the statistical significance from 
the cosmic variance efi^ect. The peak-like reionization model, 
on the other hand, has some distinct features in the shape of 
ionization fraction, and of the polarization power spectrum 
as well. The most pronounced manifestation of the peak-like 
reionization model is the localized features in the polariza- 
tion power spectrum which difi'ers from the standard single 
reionization model. We reckon that such kind of deviation 
from the standard reionization model, in case of confirma- 
tion by the Planck data, will bo significant for investigation 
of unstable particles or any relic decaying during the 'dark 
age' of the Universe. 

Note that in this paper we do not consider the sec- 
ondary anisotropics and polarization produced by the peak- 
like reionization at high redshifts. These off'ccts seem to bo 
important if we take into account the fact that the relax.- 
ation of the peculiar velocity of baryonic matter and dark 
matter at z ~ 200 is completed and we can have the analog 
of the Ostriker-Vishniak effect and the Doppler effect, but 
for specific shape of the ionization fraction. These effects will 
be investigated in the next paper. 
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ance limit for the Planck missions. 
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model 5b and the standard at Zreion — 6 and the dash-dot line 
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dotted (/sky = 0.65) and solid lines (/sky = 1) represent the 
cosmic variance limit for the Planck missions. 
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